Boson-fermion stars: exploring different configurations 
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^ ! I. INTRODUCTION 

o . 
o ■ 

■ Scalar fields have been playing an increasingly important role in recent investigations in cosmology and astrophysics. 
-~ | A number of such fields have been considered, ranging from the dilaton in superstrings, to the Higgs field and from the 
postulated inflaton, to the axion. This is one of the reasons for the interest in the properties of stellar objects made 
out of boson fields, called boson stars, since it is quite possible that such objects play, or have played, a considerable 
role as components of the dark matter of the universe. Moreover, such exotic objects provide us with unique insights 
into very high energy physics which, apart from the imprints left in the Early Universe, is currently almost completely 
\ unaccessible to experiment pj. 
J> ■ Boson stars can be seen as macroscopic boson condensates, formed by the vacuum expectation value of scalar fields, 
| the particles being the excitations over the vacuum expectation value. First studied by Ruffini and Bonnazola 0, 
their work was later generalised to include self- interactions gauge charges 0, excited solutions of the coupled 
Einstein and Klein-Gordon equations, as well as non-topological soliton solutions |a, |(J , Q-stars 0, oscillating soliton 
stars 111, and even stars with two bosons with different masses 0. The stability of these stars was studied in |Tfij 
and [ll||. Two complete reviews on the subject can be seen in references ^3 an d [13 • 

However, if these stars had their origin in some primordial gas of both bosons and fermions, one would certainly 
expect to find objects made out of a mixture of these two types of particles. Even if they were purely bosonic or 
fermionic at their origin, they would later be susceptible to a considerable amount of contamination by bosons or 
fermions, respectively. This is why it is important to investigate these mixed objects, now called boson-fermion stars. 
We shall not study here the important question of how these stars, both pure boson or boson fermion might have 
been formed, by gravitational condensation, out of a primordial gas. This issue has been addressed in part, for pure 
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scalar fields, in 14] and [15l |. We shall assume these mixed stellar objects to exist and study the configurations they 
. . , can take as stable systems. 

Their stability and main properties have been investigated in [Hi Il7l Il8| . Generalising simple heuristic arguments 
by Thirring it was then pointed out that many of the general properties could be understood on the basis of 
an interplay between the Pauli principle and the Heisenberg principle, the first being ultimately responsible for the 
stability of the fermionic component and the uncertainty principle being responsible for the stability of the bosonic 
component. The dependence on these two fundamental principles of quantum mechanics is responsible for the very 
different characteristics shown by these stars, when we have configurations that are either predominantly fermionic 
or predominantly bosonic. 

In most of the former investigations on the subject, the fermions were taken to be neutrons. We now relax this 
assumption, taking the fermion mass as a free parameter, accepting as reasonable that we may have different varieties 
of fermions in a primordial gas. As for the boson field, we take it to be a real, massive, self-interacting scalar field, 
with the mass of its associated particle also taken as a free parameter. 
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Now that new ideas, like the one provided by fermion balls |20j | and supermassive boson stars |2l| . are being put 
forward to explain the properties of galaxies, particularly the possible existence of supermassive dark objects at galactic 
centers, or of using boson stars in the role of gravitational lenses |22j it is appropriate to come back to the question 
of fermion-boson stars and to investigate if they can provide us with alternative scenarios of interest in astrophysics 
and cosmology, in particular when coupled with the problem of dark matter. As a result of the properties of bosons 
and fermions, we are able, for instance, to construct compact fermionic cores, more or less massive, depending on the 
mass taken by the fermion, surrounded by immense halos of bosons, whose sizes are also dependent on the mass of 
the corresponding boson. Very different objects can be constructed, by varying the masses of their components. 

Modern high-energy physics theories providing us with a large variety of bosons and fermions as serious candidates 
for dark matter, including massive sterile neutrinos, gravitinos, axinos as well as Higgs fields, axions, heavy scalars 
resulting from non-thermal decays of inflatons, and still others 23], make it worth while to pursue these theoretical 
constructions and explorations. 

II. THE MODEL 

The equations have been derived before, nevertheless, and for completeness, we repeat them here. The metric we 
shall use is the one appropriate to a stationary, spherically symmetric distribution of matter: 

ds 2 = -B(r) dt 2 + A{r) dr 2 + r 2 d0 2 + r 2 sm 2 9d(/) 2 ; (1) 

The Lagrangian for the massive real scalar field is defined by 

C=- l -g^d^d^- l -m 2 V 2 (2) 

and from it we derive the energy-momentum tensor: 

Tpv = d^epdvp - g^ u (J^g^dpipdaLp + ^m 2 ip 2 + jA<^ 4 ^ • (3) 

We expand the scalar field in terms of creation and annihilation operators in the usual manner |lfil llfil ITtI ] , in such 
a way that tp + = tp, as is appropriate for a real field: 

<p(r, t) = E n 1 /2 (a„^„(r) e — »' + a+^+(r)e- 4 ) (4) 

with 

[a m , ot n ] = S mn . (5) 

We shall assume in our problem all the bosons to be in the discrete quantum state defined by n = 0, defining the 
lowest energy of the star, luq. This corresponds to a spherically symmetric and nodeless wave function fo(r). It also 
means that the bosons, in number Nb, will be in the quantum state 

\N B ) = (N B \)- 1 J 2 (a+) B T |0>. (6) 

Following Breit et al. [24|, we find that the canonical commutation relations require that 

i7rr 2 ^\ip \ 2 dr=l. (7) 

The Einstein equations being classical equations, in order to write their right hand side we take the following expec- 
tation value for the energy-momentum tensor: 

V = (N B \ :T„ v : \N b ) - d^cp+d^ - ^ 9fw [gP° d p <p+ d a tp c + m 2 ip+y c + ^Mvc} 4 ^ (8) 

where we introduced, with the help of eq. 0, the semi-classical field 

<p c (r, t) = M 1/2 (^b + l) 1/2 <Po(r)e- Wot ■ (9) 



The normalisation condition now takes the form 

WI^^^ + s)- (10) 

It is well known that Fermi statistics leads the energy-momentum tensor to approach the perfect fluid form, when 
the number of fermions present is large 0, l25| . The equation we shall use for the degenerate fermion gas was first 
derived by Chandrasekhar and has the following parametric form: 

p = K{smkT-T) (11) 

K T 
p = — (sinhT-8sinh- + 3T), (12) 
3 2 

the constant K being given by K = m^/(i2i: 2 ) where my is the mass of the fermion, which we shall leave as a free 
parameter. The parameter T is defined in terms of the maximum value of the momentum of the Fermi distribution 
at radius r by 




T(r)=41og f^ + Jl + f^ . (13) 



Finally, the particle density n(r) is given by the expression 

3tt 2 3tt 2 Smn 4 

In the absence of an explicit non-gravitational interaction between fermions and bosons, the Einstein equations are 



n W ^ = -Lsinh 3 ^. (14) 



R pv - 2^^ R = 8?rG ( T ^(fc) + T^(p, p)) , (15) 

with T^(ip c ) given by eq. JHJl and the fermion part is defined by the perfect fluid expression 

T ixv{p,p) = (p + p)Ufj,u u + g^p (16) 

being the 4-velocity of the fluid. From the Bianchi identities we find that the covariant derivative of the energy- 
momentum tensor must be zero, from which condition we derive the following two equations, which will then be 
coupled to the Einstein equations: 

■dr = -2^ {P + P) (17) 

and, for the scalar field, 

Dip c — m 2 (p c + A |v? c |Ve (18) 

where □ is the d'Alcmbcrtian operator associated with our metric. Before we proceed, we must remember that the 
explicit time dependence that appears in this equation and in eq. (JSJ disappears, the exponential factor exp(— luj^t) 
being cancelled out; the result of the action of time derivatives is then the term uj 2 ip c . 

We introduce some redefinitions of the variables [17| , which enable us to write the equations in a very convenient 
dimensionless form (c = h = 1): 

X = mr; a{x) = (4.G) 1 / V(r, 0) ; SI = ^ ; A = ; p = ^ ; p=^. (19) 



m 

The equations to be integrated become: 



A' = - A(A 1] + x A 2 (tf^ + ^+a 2 + l\^ + 2p) ; (20) 
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B(A - 1) A ( a a 2 a' 2 2 1 , \ 
B = — '-+xAB(n 2 — + — -o 2 - -\o A + 2p\ ■ (21) 



where the last equation is derived from eq. I|17fl arid 



p = a (sinhT - T) , (24) 

1 / T T\ 

p = -a sinhT-8sinh- + - (25) 



and a — m^/(8nrn 2 m 2 Pl ), mpi — G x l 2 being the Planck mass. 
To calculate the total mass of the star we use 

A(*) = (l-2^)"\ (26) 

Integrating the equations up to a point sufficiently far out that we may consider the fermion and boson terms to be 
negligible, we obtain the mass from x 

*(«) = f(l-i). (27) 

The number of fermions and the number of bosons in our star can also be obtained from analytical expressions. 
From eqs. (|1()|> and l|14|l we find 

N B = lu f^xf^-ipldr^nf 1 ^) 2 \J^-a 2 x 2 dx (28) 
J V B \ m ) J V B 

J Airr 2 y/Andr = i- {^f J sinh 3 | J~Ax 2 dx , (29) 



N f 

where we have written both the original expressions and the ones obtained after the redefinitions introduced above. 
These values enter into the calculation of the binding energy, defined by 



E B = M - (m b N b + m f N f ) . (30) 

A positive value of E B means an unstable configuration. For a more complete discussion of the issue of stability in 
boson- fermion stars we refer the reader to |Tfij . 

Although we cannot separate the mass of the fermionic component of the star, we introduce in addition, as a useful 
approximation, but only as an approximation, the following quantity for the fermion core: 



M f = J 4nr 2 pdr. (31) 

From the structure of the equations, we find that the initial conditions are subjected to the following restrictions: 
A(0) = 1 and er'(0) = 0, with cr(0) and T(0) free parameters. There is no restriction on B, except that B = 1 at 
infinity. As the equations are linear in B, this can always be obtained by rescaling B. Apart from B, only Q needs 
to be rescaled. 

Using the redefined variables, once we fix the initial values ip c (0) and T(0), and the value of a, our equations and 
their solutions are mathematically the same, for any combination of the fermion and boson masses giving the same a 
in eqs. 1)24(1 and l|25|l. Only the scales of the problem will be changed. From eq. I(19|l we see that, putting <t(0) = 1, 
¥> c (0) is close to the Planck energy. 

In the numerical simulations that follow we put A = 0. The dynamical influence of this constant on the total mass 
of the star has been explored in references ^3 and but such influence is not particularly relevant for our 
considerations here. We could also have introduced direct couplings between the fermion and the boson fields, either 
in the spirit it has been done in an( i US or ! m a more phenomenological way, in (2?|]. They lead to limits in 
the couplings and masses of the scalar particles. We made sure that, using values within the limits set up in [27]], no 
important alterations would came to our results. 



III. NUMERICAL SIMULATIONS AND DISCUSSION 



From the equation for the scalar field, comparing the mass term with the term involving the eigenvalue, correspond- 
ing to the second time derivative of tp, we would in general expect the eigenvalue w to be approximately of the same 



order of magnitude as the mass of the scalar field. This is confirmed by our numerical simulations. The eigenvalue 
to appears as a factor in front of eq. (|28(l for the number of scalars. In most simulations, the mass of the scalar field 
is indeed very small. No wonder then that, to obtain configurations with a reasonable boson population, we find it 
necessary to use at the origin values for ip c (0) which are close to the Planck value. 

For the fermions the situation is different. Using a value of T(0) around 1, and the neutron mass (m n ) for the 
fermion, we find parameters typical of neutron stars, far away from the Planck scale. When we change my, these 
values are obviously scaled accordingly. 

All we have to ask now is how reasonable is it to take such large values for <^ c (0) in the simulations. If we assume 
such stars to have formed in the early universe, as usually is the case when we deal with objects like primordial black 
holes, then planckian-like values seem reasonable. When the stars are formed at later stages, this assumption is more 
questionable. It is thus important that we do not restrict our attention to a too small domain of values for ip c (Q). 
Likewise, we shall explore a large domain of masses for the fermions and the bosons, not restricting ourselves to known 
or expected masses. 

We shall apply an iterative shooting method to determine the eigenvalue to. A technical problem appears when 
we deal with a combination of masses such that the boson wave function quickly decreases to very small values, well 
inside the fermion core. We have to make sure this is not due to an overshooting. The way to do it is to make sure 
that, by displacing the node of <j(x), there are no changes in the total mass of the star and in the number of bosons, 
which in turn requires that the eigenvalue be fixed with great accuracy |l7^ . In what follows we shall study two 
different situations which differ mainly by the relative magnitude of the bosonic and fermionic terms in the Einstein 
equations. 

A. Fermionic and bosonic terms of the same order of magnitude 

To set the problem, we begin by what looks like a natural choice, taking the fermion mass equal to the neutron 
mass and assuming the different terms in our equations to be of the same order of magnitude. This requires a tp c (0) 
not very distant from the Planck scale. We will now analyse four different cases: 
Case 1: T(0) = 1 ; <p c (0) = 2.3 x I0 17 GeV (<r(0) = 0.1) ; m f = m„ ; m b = 2.71 x ICT 20 GeV 

In this case we find 

p(0) = 9.96 x 10 13 gcnT 3 ; 
M = 4.76 x I0 33 g « M Q {M f = 3.44 x I0 32 g « O.1M ) ; 
N b = 9.32 x I0 76 (N b m b = 4.50 x I0 33 g w M ) ; 
N f = 2.08xI0 56 (Nfirif = 3.47 x I0 32 g ~ O.lAf©) ; 
E B /c 2 = -9.0 x I0 31 g; 
r c = 15.9 km; 
rgo = 68.9 km; 

where p{0) is the central density corresponding to the fermionic component, r c is the radius of the fermion core and 
rgo the radius within which 90% of the total mass of the star is concentrated. This is a value we shall always quote 
whenever, as is the present case, the tail of the wave function of the scalar field extends beyond the fermion core, 
forming a halo around it. From the inspection of N b m b and Nfirif, or M and Mf, we conclude that the present 
configuration is dominated by the boson component which accounts for more than 90% of the total mass. This is an 
example of values we have found in ref. |l7j |. The binding energy was converted into the mass equivalent. 

Keeping the same fermion mass, we go to higher boson masses by giving smaller values to a. For instance, with 
a = 10 -24 , we find m b = 2.7 x I0~ 8 GeV and we have to decide what is a reasonable choice for <p c (0). One possibility 
is to assume that, once the gravitational field couples to the energy, the terms in the equations are of the same order 
of magnitude, which was the philosophy we used above, except that this time this requires a value for y c (0) « 10 6 
GeV. It is instructive to compare this with the case where we give to the scalar field a value close to the Planck scale. 
First the case where the terms are about equal. 

Case 2: T(0) = 1 ; <p c (0) = 2.3 x I0 6 GeV(cr(0) w I0" 12 ) ; m f = m n ; m h = 2.71 x I0~ 8 GeV 
Now we get: 

p(0) = 9.96 x I0 13 gcm- 3 ; 
M = 6 X I0 32 g re O.1M (M f « M) ; 
N b = 3.70 x I0 47 (N b m b = 1.79 x 10 16 gw I0" 17 M Q ) ; 



N f = 3.62 x 10 56 {N f m f = 6.04 x 10 32 g w 0.1M©) ; 
£ B /c 2 = -5.44 x 10 30 g; 
r c = 20.8 km; 

where rgo is not quoted, as now the scalar wave function goes practically to zero well inside the fermion radius. The 
star is entirely fermion dominated and it has the corresponding typical physical quantities, like mass and radius. If 
we increase T(0), the number of fermions increases, as expected, but the number of bosons remains practically the 
same. This is a typical configuration to be considered when dealing with neutron stars contaminated by WIMPS. If 
the masses of the WIMPS are of the order of 10 -5 eV or above, they will concentrate in the core of the star. 

In the following example, the scalar field is closer to the Planck scale: 
Case 3: T(0) = 1 ; <p c (Q) = 2.3 x 10 17 GeV (cr(0) = 0.1) ; m f = m n ; m b = 2.71 x 10~ 8 GeV 

Now we find 

p(0) = 9.96 x 10 13 gcm" 3 ; 
M = 5.23 x 10 21 g = 10~ 12 M Q (Mf = 3.14 x 10" 3 g) ; 
N b = 1.11 xlO 53 (N b m b = 5.33 x 10 21 g « 10" 12 M Q ) ; 
N f = 1.89 xlO 21 {Nfirif = 3.17 x 10" 3 g); 
E B /c 2 = -9.99 x 10 19 g; 
r c = 3.13 x 10~ 6 cm; 
r 90 = 7.0 x 10~ 6 cm . 

We have a configuration which is entirely different from the preceding case, although we have kept the same central 
density p(0). It is now a typical soft boson star, with a microscopic radius, but a still respectable mass. The fermion 
component has almost completely disappeared. The increase in <p c (0) has resulted, not only in a substantial increase 
in Nf,, but also in an enormous reduction in the number of fermions. No wonder this has happened as, with the choice 
made for ip c (0), the terms in the equations associated with the bosons are very much larger than the terms associated 
with the fermions. In this sense, it is not a natural choice. 

If we increase T(0) we find that, around T(0) = 3, the configuration will suddenly become fermion dominated, a 
property that has been noticed and investigated in the work of reference 01 1 an d we shall have a configuration similar 
to the one in case 2 above, albeit with a very different value for C (O). 

Before going to a different kind of limit, we mention one case corresponding to the fermions and bosons having 
much higher masses than the ones used so far: 

Case 4: T(0) = 1 ; ip c (0) = 2.3 x 10 6 GeV O(0) » 10~ 12 ) ; m f = 10 7 GeV ; m b = 3.08 x 10 6 GeV 
In this case: 

p(0) = 1.29 x 10 42 gcm- 3 ; 
M = 5.28 x 10 18 g « 10~ 15 M Q (M f « M) ; 
N b = 2.90 x 10 19 {N b m b = 1.59 x 10 2 g) ; 
Nf = 2.99 xlO 35 (Nfmf = 5.33 x 10 18 g) ; 
E B /c 2 = -4.79 x 10 16 g; 
r c = 1.83 x 10 _8 cm 

It is really an object of atomic size, with characteristics close to primordial black holes, without being one. The value 
we took for <p c (0) is within the same order of magnitude as the masses and is a case where the fermion and boson 
terms in the equations are of the same order of magnitude. 

Again, it is completely fermion dominated. Except for the case 1 above, in all the others, whenever the terms 
associated with the bosons and the fermions are of about the same order of magnitude, the resulting configuration is 
dominated by fermions. This happens, for rrif — m n , when m b is larger than 10 -5 eV, as in case 2 above. 

B. Fermionic terms dominate 

We shall now investigate what happens when m/ is equal to or smaller than the nucleon mass and the boson mass is 
really very small. This means that we have to increase a. The fermion terms will necessarily dominate the equations, 
whatever value we use for <p c {0), if we keep it, as we shall do, equal or below the Planck scale. However, due to the 



smallness of the boson mass, this does not mean that the star will be fermion dominated. In case of pure boson stars, 
we know that the total mass of the star tends to change according to 1/wi 2 - As a first example we take: 
Case 5: T(0) = 1 ; <p c (0) = 2.3 x 10 17 GeV (cr(0) = 0.1) ; m f = m n ; m b = 8.56 x l(T 23 GeV 
Now we have 

p(0) = 9.96 x 10 13 gcm~ 3 ; 
M = 1.65 x 10 36 g « 10 3 M e (M f = 5.99 x 10 32 g » 0.1M Q ) ; 
N b = 1.10 x 10 82 (N b m b = 1.68 x 10 36 g » 1O 3 M ) ; 
A// = 3.61 x 10 56 (Nfinf = 6.04 x 10 32 g « 0.1M®) ; 
£ B /c 2 = -3.16 x 10 34 g; 
r c = 20.8 km; 
r 90 = 2.18 x 10 4 km . 

We have a complicated configuration, with a core typical of a neutron star, and a very massive boson halo extending 
well beyond. The bosons dominate the total mass by many orders of magnitude. 

Had we used a still smaller boson mass, for instance m b — 2.71 x 10~ 32 GeV, we would have got a star with a total 
mass of M = 5.26 x 10 45 g s» 10 12 M Q and a halo extending up to 6.92 x 10 13 km = 2.24 pc, most of the mass coming 
from the boson component (the approximate value for the fermion component is Mf = 5.99 x 10 32 g w 0.5Mq and 
r c = 20.8 km, the same as above). 

The authors of [2fJ have tried to show that, choosing a fermion mass of about 16 KeV (for instance, a sterile 
neutrino), it was possible to explain the main dynamical features of the supermassive compact dark objects, known 
to exist in the Galactic center and in the center of M87, with the help of what they have dubbed fermion balls, 
which they propose as an alternative explanation to the more usual black hole model. Better observations have, 
however, ruled out this model, but still leaving opened the possibility of a supermassive boson star plj . 

These ideas suggest that we consider a boson-fermion star with both a small fermion mass and a very small boson 
mass. When we consider this possibility we get the following configuration: 

Case 6: T(0) = 1 ; <p c (0) = 2.3 x 10 17 GeV (cr(0) « 0.1) ; m f = 10" 5 GeV ; m b = 3.08 x 10 - 30 GeV 
For this last case we have 

gcm~ 3 ; 

4.18 x 10 43 g w 10 10 Af o (M f = 3.03 x 10 42 g w 10 9 M Q ) ; 
(N b m b = 3.96 x 10 43 g w 10 10 M o ) ; 
(Nftrif = 3.05 x 10 42 g w 1O 9 M ) ; 

L g; 

cm = 4.53 x 10~ 3 pc 
cm = 1.97 x 10~ 2 pc 



P(0) 


= 1.29 x 10~ 6 


M 


= 4.18 X 10 43 ; 


N b 


= 7.22 x 10 96 


N f 


= 1.71 x 10 71 


E B /c 2 


= -7.98 x 10 4 


r c 


= 1.40 x io 11 : 


rgo 


= 6.09 x io 11 : 



The fermion core and the boson halo form a supermassive compact object with the halo extending up to 2 x 10~ 2 
pc. Keeping the same fermion mass, but taking, for instance, a boson mass in the range of 10~ 32 GeV, we would get 
a total mass of the order of 10 13 solar masses, with a halo extending up to 6 pc, while the fermion component would 
have practically the same mass and the same radius as in case 6. 

Would it in principle be possible to distinguish a boson-fermion system from a fermion ball? Once again this could 
be done through the careful investigation of the orbits of stars around a candidate object. We must take into account 
that the mass density profile of a boson-fermion configuration, with an extended boson halo, would be different from 
the profile of a fermion ball. This should be reflected in the dynamics of the orbits of stars inside such a halo, exactly 
as there appears a difference between the fermion ball and the black hole models 20] . We are currently investigating 
this possibility and the results will be reported in a separate publication [2flj . 

We are thus in the presence of a system which can take many different configurations, characterised by very different 
masses and radii, ranging from objects of atomic size and masses of the order of 10 18 g, to objects having galactic 
masses and extending up to a few light years. To show this flexibility was the main purpose of the present work. 
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